First direct measurement of optical phonons in 2D plasma crystals 
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Spectra of phonons with out-of-plane polarization were studied experimentally in a 2D plasma 
crystal. The dispersion relation was directly measured for the first time using a novel method of 
particle imaging. The out-of-plane mode was proven to have negative optical dispersion, comparison 
with theory showed good agreement. The effect of the plasma wakes on the dispersion relation is 
briefly discussed. 
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Complex plasmas consist of particles immersed in a 
weakly ionized gas-discharge plasmas [H, [3, Due to 
the absorption of ambient electrons and ions, micropar- 
ticles acquire significant (negative) charges and can form 
strongly coupled systems analogous to conventional soft 
matter (complex fluids). Two-dimensional (2D) com- 
plex plasmas are particularly convenient systems for the 
detailed experimental studies of strongly coupled phe- 
nomena at the atomistic level. The essential advan- 
tage of such systems is that one can gain complete in- 
formation about the state of the entire particle ensem- 
ble in the phase space, which is an invaluable advantage 
for the investigation of collective processes occurring in 
strongly coupled media, such as transport, phase transi- 
tions, structural quenching, etc [H, 

2D complex plasmas are usually studied in ground- 
based experiments with RF discharges. Microparticles 
injected in such discharges levitate in the sheath region 
near the bottom electrode, where the electric field can 
balance gravity. Under certain conditions the particles 
can form a monolayer, and if the electrostatic coupling 
between them is strong enough, particles arrange them- 
selves into ordered structures - 2D plasma crystals [1] . In 
such crystals, numerous collective phenomena like melt- 
ing and recrystallization, mass and heat diffusion, soli- 
tons and shocks SBSBSIii, etc. have been reported. 
Propagation of dust-lattice (DL) waves in 2D lattices is 
often used as a diagnostic tool to determine parameters 
of the plasma crystal |Tl[ . 

In 2D complex plasmas, like in any (strongly cou- 
pled) 2D system, two in-plane wave modes can be sus- 
tained (here we naturally leave aside polydisperse mix- 
tures). In crystals, both modes have an acoustic dis- 
persion, one of them is longitudinal, another is trans- 
verse. Since the strength of the vertical confinement 
in such systems is finite, there is a third fundamental 
wave mode associated with the out-of-plane oscillations 
[13, [H 0, [H, [B, [13, [13. Theory predicts that this 
mode has a negative (or inverse) optical dispersion and 
depends critically on the parameters of the plasma wakes 
(that are formed downstream from each particle levitat- 
ing in the sheath [3, [3]). So far, the systematic ex- 
perimental studies of the DL waves were limited by the 
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FIG. 1: (Color online) Sketch of the experimental setup. Mi- 
croparticles are confined above the RF electrode and are illu- 
minated with a horizontal laser sheet having a Gaussian pro- 
file in the vertical direction. The monolayer is levitated well 
below the peak of the laser intensity, which results in strong 
intensity variations of the scattered light upon the vertical 
displacement of individual particles. The upward/downward 
displacement corresponds to positive/negative intensity varia- 
tion. Inset shows elementary cell of the hexagonal lattice and 
the frame of reference chosen in this paper, the orientation of 
the wave vector k is measured in respect to the x axis. 



in-plane modes, whereas direct measurements of the out- 
of-plane dispersion relation have never been carried out 
(see, however, Refs. [l^, [HI). To a large extent, this 
is due to the lack of reliable diagnostics of the vertical 
dynamics of individual particles. 

In this Letter, we report on the first direct measure- 
ments of the out-of-plane DL wave mode in a 2D plasma 
crystal. The dispersion relation is recovered in a broad 
range of wave vectors, revealing remarkably good agree- 
ment with the existing theory. The measurements were 
performed by employing a new detection technique of 
the vertical (out-of-plane) particle displacement. The ob- 
tained results provide us with new insights into the mech- 
anisms governing the individual particle dynamics in 2D 
complex plasmas (in particular, highlight the important 
role of the plasma wakes) and allow us to specify the 
conditions when such systems can be used to investigate 
generic properties of (classical) 2D solids and liquids. 

A sketch of the experimental setup is shown in Fig. [T] 
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We used a capacitively coupled RF glow discharge at 
13.56 MHz. The argon pressure p was between 0.5 Pa 
and 1 Pa and the RF peak-to-peak voltage Vpp was be- 
tween 175 V and 310 V (which corresponds to a forward 
RF power P between 5 W and 20 W) . The plasma param- 
eters in the bulk discharge were deduced from Langmuir 
probe measurements, yielding the electron temperature 
Te = 2.5 eV and the electron density rig = 2 x 10^ cm~'^ 
at p =0.66 Pa and P=20 W ji^l. A 2D particle sus- 
pension was formed by levitating melamine formaldehyde 
microspheres in the sheath above the RF electrode. The 
particles had a diameter of 8.77 ± 0.14 /im and a mass 
m = 5.3 X 10~^^ kg. The diameter of the obtained crys- 
talline structure was about 50-60 mm, depending on the 
number of injected particles. The microparticles were il- 
luminated by a horizontal laser sheet which had a Gaus- 
sian profile in the vertical direction with a standard de- 
viation cr ~ 75 /im (corresponding to a full width at half 
maximum of 175 fim). The sheet thickness was approxi- 
matively constant across the whole crystal. The particles 
were imaged through a window at the top of the chamber 
by a Photron FASTCAM 1024 PCI camera at a speed of 
250 frames per second. The horizontal coordinates x and 
y as well as velocity components Vx and Vy of individual 
particles were then extracted with sub-pixel resolution 
in each frame by using a standard particle tracking tech- 
nique (2l| . An additional side-view camera was used to 
verify that our experiments were carried out with a single 
layer of particles. 

In order to extract the vertical position z and velocity 
Vz of individual particles, we needed to employ a very 
different technique. Unlike relatively compact Coulomb 
clusters, where these values can be directly measured by 
using the side view, the particles forming large monolay- 
ers cannot be individually resolved from the side. There- 
fore, the out-of-plane particle tracking can only be per- 
formed by using the top view, and then the vertical dis- 
placement should be deduced from the relative variation 
of the scattered light intensity. 

The "conventional" method of particle visualization, 
when the peak of the Gaussian intensity in the illumi- 
nating laser sheet practically coincides with the particle 
monolayer (this facilitates the in-plane tracking, espe- 
cially at higher frame rates), is not appropriate to ob- 
tain z and Vz- The reason is twofold: The first problem 
of such configuration is that the variation in the (par- 
ticle image) intensity between consecutive frames gives 
us only the magnitude of the vertical displacement (and 
hence velocity), whereas the direction remains undefined. 
Indeed, the intensity of the individual particle image at 
each frame is proportional to the local illumination in- 
tensity, which scales as 

/(z)oce-(--^— (1) 

where Zmax is the position of the intensity maximum. 
Therefore, if the mean levitation height zi^v is very close 
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FIG. 2; (Color online) Fluctuation spectra of the naturally 
excited out-of-plane DL waves in a 2D crystal. The fluctu- 
ation intensity is color-coded linearly, from dark-blue (10"'^ 
a.u.) to dark- red (1 a.u.), revealing the optical dispersion of 
the mode. The results are for a pressure of 0.8 Pa and a 
RF power of 15 W. The upper and lower panels depict the 
fluctuation spectra with the wave vector k pointed along the 
principal lattice axes, 9 = 0° and 6 = 30°, respectively (see 
Fig. [T|. The wave vector is normalized by the inverse lattice 
constant . 

to Zmax, then the magnitude of the particle displacement 
\Sz\ can be substantially larger than | 1. Ob- 

viously, a given variation of the image intensity in this 
case may equally imply either positive or negative dis- 
placement |13|. The second problem of the "conven- 
tional" method is associated with the sensitivity. Since 
Sz ^ l^icv — 2;max|, the resulting intensity variations scale 
as 6I/I ~ ^{6z/a)^. The magnitude of vertical displace- 
ment (induced by thermal fluctuations) can be estimated 
as \6z\ ~ y^Td/rridill 10 A^ni (here Td is the kinetic 
temperature of particles and Qv = 27r/v is the resonance 
frequency in their vertical confinement). This value is 
much smaller than the standard deviation cr « 75 //m. 
Thus, we conclude that the relative intensity variations 
are at the level of ^ 1 %, which makes them very difficult 
to detect. 

In our experiments we employed an alternative method 
of particle visualization. The position of the intensity 
maximum was set about 100 /im above the levitation 
height, so that |ziov — Zmax|/|'52;| w 10. This allowed us 
to overcome both drawbacks mentioned above: (i) par- 
ticles moved in the region where the vertical gradient of 
the illumination intensity remained positive (which al- 
lowed us to avoid ambiguity in determining the direc- 
tion of the displacement), and (ii) the relative intensity 
variations scaled almost linearly with the displacement, 
61 /I K |ziev — -^maxl^^/c^, SO that the resulting magni- 
tude oi SI/I was about 15 times larger than that in the 
"conventional" tracking configuration. It is worth noting 
that this alternative visualization method does not crit- 
ically affect the quality of the tracking in the horizontal 
xy plane compared to the "conventional" tracking con- 
figuration. A separate test has indeed shown that the 
pixel locking effect is roughly at the same level and that 
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FIG. 3: (Color online) Comparison of the experimentally mea- 
sured DL dispersion relations with theory. Shown are the 2D 
in-plane and out-of-plane dispersion relations deduced from 
the experiments (at RF power of 15 W and pressure of 0.8 Pa) 
and the theoretical curves. The results are presented for two 
directions of the wave vector k (for 6 — 0° and 6 — 30°, see 
Fig. [ij . The dashed lines are theoretical dispersion relations 
for pure Yukawa interaction between particles (with the same 
parameters as in the experiment, see Table |T]), the solid lines 
represent the case when the interaction with plasma wakes is 
taken into account. The shown range of |k| is limited by the 
first Brillouin zone and therefore depends on the direction: 
|k|a = |7r for 6 = 0° and |k|a = ^vr for 61 = 30°. The wave 

vector is normalized by the inverse lattice constant a^^ . 



the in-plane fluctuation spectra are similar. 

By performing the Fourier transform in time and space 
of the determined particle velocities, the dispersion rela- 
tions /(k) (both in-plane and out-of-plane) can then be 
obtained for a 2D crystal. 

The resulting fluctuation spectra of the out-of-plane 
waves are shown in Fig. [3 for two principal orientations 
of the wave vector, at 6' = 0° and 30°. The spectra 
represent the wave energy distribution in the (/, k) space, 
so that the "ridge" of this distribution yields the wave 
dispersion relation. The most conspicuous feature of the 
measurements is the optical character of the dispersion 
relation - the frequency of the long-wavelength waves is 
finite. Moreover, the dispersion at long wavelengths is 
negative, i.e., the wave frequency falls off as the wave 
number increases. At larger |k|, the wave dispersion is 
different for 6* = 0° and 30°. 

Experimental observation of the out-of-plane wave dis- 
persion is the main result of our paper. Figure [2] shows 
that the wave frequency as the function of the wave vec- 
tor k changes in a narrow interval of 17 Hz 20 Hz, 
which makes resolving the wave dispersion a challenging 
task. Our method, however, allowed us to clearly identify 
the wave dispersion. Indeed, the frequency resolution of 
sa 1.5 Hz (defined as the standard deviation of a Gaussian 
fit) is smaller than the measured wave frequency range 
A/(« 3 Hz). 

Our method allowed us to measure all three wave 



TABLE I: Parameters of a 2D crystal extracted from different 
experiments (marked from I to IV). The lattice constant a 
is obtained from the position of the first peak of the pair 
correlation function. Lattice parameters Zd and k = a/\ are 
deduced from the in-plane phonon spectra with accuracy « 
15 % and ~ 30 %, respectively. The accuracy for the measured 
vertical resonance frequency /v is !5i 8 %. The frequency range 
A/ is the difference between /v and the minimum frequency 
of the out-of-plane mode at 6' = 0° (see Fig. [TJ . Subscripts 
"exp" and "th" stand for "experimental" and "theoretical", 
respectively, A/jj^ is calculated using the theory of Ref. ;18| . 
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0.6 


0.4 
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16.3 
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3.9 



modes of a 2D plasma crystal in a single experimental 
run. In Fig. [3] we show the fluctuation spectra of out-of- 
plane and in-plane wave modes obtained from the same 
experiment as in Fig. [2l for two principal lattice orien- 
tations, 9 — 0° and 30°. Only the flrst Brillouin zone is 
shown here. Lines represent theoretical dispersion rela- 
tions which depend on two lattice parameters - particle 
charge number Zd and the lattice screening parameter 
K — a/ A, where a is the interparticle distance (lattice 
constant) and A the effective screening length. The theo- 
retical dispersion relations were calculated using the fol- 
lowing procedure. First, the lattice parameters were de- 
duced from the low-|k| part of in-plane fluctuation spec- 
tra using the method of Ref. The obtained values 
of Zd and k are shown in Table H] for four different ex- 
periments. Second, these values were used in the theory 
of Ref. fl§| to calculate the whole dispersion curves (the 
difference between the dashed and solid curves will be 
discussed later). 

As can be seen in Fig. [3l all three branches of the ex- 
perimental wave spectra can be simultaneously fltted by 
theoretical curves calculated with the same set of param- 
eters. This provides a useful consistency check for our 
experimental method and also for theory. Good over- 
all a gree ment of our measurements and the theory of 
Ref. [18| serves as a verification that our method cor- 
rectly measures the out-of-plane motion of particles and 
also as an indication that the theory of Ref. 3] correctly 
captures the essential physics of 2D plasma crystals. 

The out-of-plane wave mode can also be useful for diag- 
nostic purpose. The lattice parameters Zd and k are usu- 
ally determined by fitting in-plane dispersion relations 
Another approach is to measure the propagation 
speed of externally excited pulses [lH . But under certain 
conditions (e.g., near the melting line) it is not practical 
and one must rely on natural waves and, in this case, 
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in-plane fluctuation spectra can be very noisy. There- 
fore, the additional fitting of the out-of-plane dispersion 
relation provides an important complementary tool to 
check consistency of the measurements. As an exam- 
ple, in Table U we compare the experimental frequency 
range for the out-of-plane mode A/oxp to the theoreti- 



cally calculated one A/th [18|. One can see that exper- 



imental and theoretical values agree well in most cases 
(except of experiment III that we discuss later), which 
is a strong indication that the crystal parameters were 
measured properly. 

Next, we discuss how the out-of-plane waves are af- 
fected by the size of the plasma crystal. Theory 3, 15, 
18, suggests that the frequency of the out-of-plane 
mode varies between /v and /v — A/, in a narrow fre- 
quency range A/. It scales with the lattice constant as 
A/th C!C r2|)L oc a~^, where rip^ = e'^ Z^/mda^ is the 
(squared) frequency scale of DL waves (it is assumed 
that 27r/v ^ J^dl)- The exact factor that determines 
the magnitude of A/th has a relatively weak dependence 
on K (which, in turn, only varies within 30 % in our ex- 
periments) and can be obtained from the rigorous theory 
for a 2D crystal [l^. Thus, for bigger crystals where the 
lattice constant a is smaller (experiments I, II, and IV in 
Table H]), the pho non frequency range A/ is larger, which 
makes it easier to resolve in an experiment. In addition, 
in big crystals the central part that is used for data anal- 
ysis is fairly homogeneous. In contrast, for small crystals 
(experiment III in Table [l| the situation is the opposite. 
The frequency range A/ oc a^^ is small and the fluctua- 
tion spectrum is broadened due to the significant varia- 
tion of a within the analysis area. This clearly results in 
greater discrepancy between experiment III and theory, 
up to 50 %. Thus, the out-of-plane wave mode can be 
better studied using bigger crystals. 

It is worth mentioning that theory 3, Hm predicts 
that both the out-of-plane and in-plane modes are af- 
fected by the presence of plasma wakes. If we introduce 
the effective positive charge q accumulated in the wake 
(normalized by the absolute value of particle charge), 
then the frequency of the in-plane modes is diminished by 
a factor ^/l — q, whereas the frequency range of the out- 
of-plane mode. A/, is proportional to (1 — q). Figure [3] 
illustrates this effect: The dashed curves are calculated 
for a crystal with a pure Yukawa interaction between par- 
ticles, whereas the solid curves are for the case when the 
additional particle-wake interactions are taken into ac- 
count (with q = 0.2) [ll[ll[ll|. This "wake correction" 
apparently improves agreement with the experiments. 

Direct observations of the atomistic dynamics in "reg- 
ular" 2D or quasi-2D systems (such as carbon nanotubes 
[i^ l and graphite are inhibited, and therefore ex- 

perimental measurements of phonon spectra in such sys- 
tems rely on indirect methods, e.g., Raman scattering 
technique 23|, |2j|. Natural model systems, such as crys- 



talline complex plasmas, provide us with a unique op- 
portunity to study generic properties of wave modes sus- 
tained in classical strongly coupled ensembles, explore 
limits of small "nano" -clusters, investigate peculiarities 
of the atomistic dynamics upon nonlinear mode coupling, 
etc. As compared to previously reported investigations of 
wave modes in 2D plasma crystals, the results presented 
in this Letter show a substantial improvement of the par- 
ticle imaging technique. It allowed us for the fist time to 
directly observe the optical phonon spectrum which has 
been predicted a long time ago. This development gives 
us strong grounds to believe that the direct observation 
of optical modes will allow much better discrimination of 
various subtle mechanisms affecting wave modes in 2D 
crystalline systems. 

The authors thank U. Konopka, B. Klumov, and D. 
Samsonov for valuable discussions. 
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